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Quantum gates in topological quantum computation are performed by braiding non-Abelian 
anyons. These braiding processes can presumably be performed with very low error rates. However, 
to make a topological quantum computation architecture truly scalable, even rare errors need to be 
corrected. Error correction for non-Abelian anyons is complicated by the fact that it needs to be 
performed on a continuous basis and further errors may occur while we are correcting existing ones. 
Here, we provide the first study of this problem and prove its feasibility, establishing non-Abelian 
anyons as a viable platform for scalable quantum computation. We thereby focus on Ising anyons 
as the most prominent example of non-Abelian anyons and show that for these a finite error rate 
can indeed be corrected continuously. There is a threshold error rate pc > 0 such that for all error 
rates p < pc the probability of a logical error per time-step can be made exponentially small in the 
distance of a logical qubit. 


I. INTRODUCTION 

Besides revealing spectacular features of quantum 
physics, non-Abelian anyons are sought for their poten¬ 
tial application in topological quantum computing m- 
Ising anyons are the most well-studied non-Abelian anyon 
model, since they describe the exchange statistics of lo¬ 
calized Majorana fermions [1] and are expected as ele¬ 
mentary excitations of the = | fractional quantum Hall 
state mm- A variety of condensed-matter systems have 
been proposed as potential hosts of Majorana fermions, 
see Refs. [HHII] for reviews. Ising anyons also appear 
as excitations [IMl] or ends of defect lines [THHI^ in 
several spin-lattice models. 

The set of quantum gates that can be performed topo¬ 
logically, i.e. by braiding anyons, depends on how qubits 
are encoded into the fusion space of a number of Ising 
anyons. However, the gate set will not allow for univer¬ 
sal quantum computation for any encoding [20]. In or¬ 
der to perform a universal quantum computation by use 
of Ising anyons, some gates need to be performed in a 
non-topological way pTI - [M| . Assuming that all topologi¬ 
cal operations are error-free, these non-topological oper¬ 
ations have a very high error threshold |5T] . 

Here, we want to focus on the assumption of error-free 
topological operations. It is often said that topological 
operations are “inherently fault-tolerant”. However, even 
a mass gap which is significantly higher than tempera¬ 
ture will still lead to a finite density of accidental quasi¬ 
particle excitations, and these need to be corrected if a 
scalable quantum computation architecture is to be built. 
The field of error correction for non-Abelian anyons is 
still relatively young. Error correction algorithms for 
Ising anyons |25j and other non-Abelian anyon models 
[13 mi have been benchmarked using Monte Carlo simu¬ 
lations. Ref. [H] demonstrated that even error correction 
for Fibonacci anyons can be simulated on a classical com¬ 
puter, despite them being universal for quantum compu¬ 
tation. Ref. [H] provides a threshold proof for arbitrary 
anyon models and a wide class of decoders. 

These references assume that we are able to detect all 


anyonic charges at some time, and then are able to fuse 
as many anyons as we like, without any further errors 
occurring. This picture, however, is highly idealized. In 
reality, further errors may occur while we are correcting 
existing ones. This is both due to the finite time needed 
to move existing anyons towards each other, and due to 
the fact that we will in general need more than one round 
of fusion to get rid of all existing non-Abelian anyons. 


For Abelian anyons models (such as the toric code 
[50]). it is possible to record all measurements for some 
time and only correct a net error at the final time step. 
This is a possibility we do not have with non-Abelian 
anyons: they need to be corrected “on the fly”, as the 
results of their non-Abelian braiding would be impossi¬ 
ble to unwind later on. Ref. m recently pointed out 
that even performing error correction after the comple¬ 
tion of each non-Abelian braid is not sufficient, since the 
braiding procedure will turn local errors into non-local 
ones. 


This article thus investigates the thus far unexplored 
problem of continuous error correction for non-Abelian 
models, where we focus on Ising anyons as the most 
prominent example of non-Abelian anyon. We restrict 
our study to the most trivial topological gate, the iden¬ 
tity - i.e., on the task of preserving a topologically stored 
quantum state. It is generally assumed that the thresh¬ 
olds for quantum information processing are identical to 
those for quantum information storage. Our main result 
is that a sufficiently low rate of errors can indeed be cor¬ 
rected continuously, allowing in principle to preserve a 
topologically stored quantum state indefinitely in a suf¬ 
ficiently large system. 

Sec. |H| discusses continuous error correction for Ising 
anyons and states our main result, the proof of which can 
be found in Sec. IHII 
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II. CONTINUOUS ERROR CORRECTION FOR 
ISING ANYONS 

Following Ref. [15] , we consider a square lattice of size 
L X L with periodic boundary conditions which hosts 
Ising anyons. Anyons are quasi-particles that live in 2D 
space, and this space needs to be discretized to allow for 
measurements of the local anyonic charge. Our setting 
is thus the most generic choice that is agnostic with re¬ 
spect to the physical realization of the Ising anyons. The 
chosen setting is also of interest from a coding-theoretic 
perspective, in that it provides an example of correct¬ 
ing a code without an explicit Pauli-matrix tensor prod¬ 
uct structure [25]. For proposals to realize Majorana 
fermions in nanowire hybrid systems, which are exper¬ 
imentally most advanced other geometries may 

be more adequate [SIIIID]. 

Each cell of the square lattice can carry either vac¬ 
uum 1, a, fermion if, or a (non-Abelian) anyon a. These 
satisfy the fusion rules 

Ip X if = 1 , Ip X a = a , a X a = 1 Ip . (1) 

All L X L charges are measured periodically at times 
0,1,2,.... We assume that these measurements can be 
performed flawlessly. This assumption means that we 
neglect the inevitable errors that arise during the com¬ 
plex measurements of non-Abelian anyons. However, it 
is common to study this ideal case in quantum error cor¬ 
rection to establish a foundation for later studies that 
include measurement noise. 

We study the question of whether it is possible to pre¬ 
serve a quantum state stored in this system despite a 
constant rate of errors affecting it. More specifically, we 
consider whether it is possible to preserve a certain state 
of the degenerate vacuum of the system. Transitions be¬ 
tween different vacuum states can be induced by dragging 
fermions or anyons around the torus. More realistically, 
one would consider storing a quantum state in the fusion 
space of a set of well-separated anyons |55|. However, we 
focus on the task of preserving a certain vacuum state of 
the torus for simplicity, since at sufficiently large length- 
scales, the question of correctability is independent of the 
particular encoding scheme chosen. 

Let us assume that between any two rounds of charge 
measurement a pair of fermions and a pair of anyons are 
created on each pair of adjacent cells of the lattice with 
probability p each (2p < 1). Error events can thus be 
associated with edges of the square lattice. Note that the 
case of both a pair of anyons and a pair of fermions being 
created on the same edge is indistinguishable from only 
a pair of anyons being created. We thus restrain from 
considering this case explicitly. Hopping and braiding of 
anyons can emerge from pair creation processes. While in 
reality different processes will have different rates, we can 
think of p as the probability associated with the highest- 
rate process. Finally, we assume that we are able to move 
an anyon or a fermion to an adjacent cell over the course 
of one measurement period. 


Our goal is to show that for sufficiently low rates p, 
it is possible to perform error correction (for any num¬ 
ber of time-steps) such that the probability of failure per 
error correction period is exponentially small in L. Con¬ 
sequently, a given state of the anyonic vacuum can be 
preserved for a time that is exponentially long in L. 

The basic idea behind our error correction approach 
is that it is always possible to first fuse all a anyons in 
pairs, and then fuse all ip fermions in pairs in order to 
obtain a vacuum state [25]. Note that according to the 
second of the fusion rules in Eq. Q , if a fermion and an 
anyon move to the same cell, they will fuse to an anyon. 
If that anyon is further moved around, it will carry with 
it the additional fermionic charge. We will refer to this 
informally as the anyon “swallowing” the fermion. The 
fermion will be recovered if the anyon is brought to fusion 
with an other anyon, as the fermionic parity is conserved. 
We will thus continuously fuse the anyons in pairs, recov¬ 
ering any “swallowed” fermions. Since the fermions are 
Abelian, error correction for them can be postponted: 
it suffices to correct a net error after some time. This 
is similar to the idea of updating a Pauli frame instead 
of performing actual corrections in a surface code (see 
e.g. Ref. [H]). If we are moving two anyons towards each 
other in an attempt to fuse them, further errors may hap¬ 
pen along their path that make one of them disappear. 
We will shortly discuss how we deal with this. 

In order to formally discuss error correction, we con¬ 
sider a 2 -|- 1-dimensional cubic lattice, in which time 
flows “upwards”, and with periodic boundary conditions 
in horizontal (spatial) direction. Charge measurements 
correspond to horizontal faces. It will prove convenient 
to identify error events and paths along which we move 
anyons with edges of the dual lattice of this cubic lat¬ 
tice. Error events happen between consecutive rounds of 
charge measurement and affect two adjacent cells. They 
can thus naturally be identified with vertical faces of 
the primal lattice, and, in turn, horizontal edges of the 
dual lattice. We will call the one-cell-per-time-step paths 
along which we move the anyons during error correction 
their world-lines. Error events will not be counted as part 
of this world-line. Horizontal faces of the primal lattice 
(charge measurements) are naturally identified with ver¬ 
tical edges of the dual lattice. For every charge measure¬ 
ment which detects an anyon, we consider the associated 
vertical edge of the dual lattice to be part of the anyon’s 
world-line. If an anyon is moved to an adjacent cell, 
the horizontal edge connecting the old and new vertical 
edge is also considered to be part of the anyon’s world¬ 
line. Note that error events always correspond to hori¬ 
zontal edges of the dual lattice, while anyon world-lines 
include both horizontal (intentional movements) and ver¬ 
tical (charge measurements) edges. Let W denote the 
union of all world-lines. 

Let A denote the set of all anyonic (as opposed to 
fermionic) error events and F the set of all fermionic 
error events. Both of these are subsets of the horizontal 
edges of the dual of the 2 -I- 1-dimensional cubic lattice. 
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Let dA denote the set of cells of the cubic lattice which 
have an odd number of elements of A incident upon them. 
Elements in dA correspond to unexpected changes in the 
anyonic charge, i.e., those which are not due to us inten¬ 
tionally moving an anyon to an adjacent cell. The sets 
W and dA are known to us with certainty. 

Fig. [2 summarizes the natural language definition and 
the geometrical interpretation of all symbols that are rel¬ 
evant for our proof. 

Let At and dAt denote the subsets of A and dA, respec¬ 
tively, that happen between charge measurements t — 1 
and t. Since the fusion rules in Eq. Q preserve the parity 
of the number of anyons that exist at any given time, the 
sets dAt always have even cardinality. It is the task of 
a classical error correction algorithm to form a hypothe¬ 
sis about the set At that is compatible with the given set 
dAt ■ This problem is exactly identical to the well-studied 
problem of finding the most likely error set in a toric code 
with bit-flip errors and perfect syndrome measurements 
m- We can thus employ the standard algorithm used 
to find such a pairing in the toric code case, namely an 
efficient minimum-weight perfect matching (MWPM) al¬ 
gorithm |42l |43j . The weight of a path connecting two 
lattice cells is thereby given by the 2-dimensional Li- 
norm, i.e., the Manhattan distance in the L x L lattice 
with periodic boundary conditions. We stress that de¬ 
spite the fact that we deal with anyon world-lines in 2 -|-1 
dimensions, the algorithmic part of the error correction 
problem for the anyons is a 2-dimensional one. This is in 
contrast to the error correction problem for the fermions. 

Let us call a subset of the edges of the dual lattice a 
string if there are exactly two cells of the cubic lattice 
which have exactly one of the edges incident upon them, 
and all other cells have either zero or two edges incident 
upon them. The MWPM algorithm will return strings 
that connect the elements in dAt in pairs. The union 
of these strings, which we call Ht, forms our hypothesis 
about what anyonic errors have happened between charge 
measurements t — 1 and t. Let H = [Jt^t denote the 
union of all hypotheses up to the present time. We note 
that At can in general not be decomposed into strings 
ending at elements of dAt - it can contain loops. These 
are defined as sets of edges of the dual lattice such that 
each cell of the cubic lattice has zero or two edges of the 
set incident upon them. 

Each element in dA is connected by a string which 
is a subset of iL to a contemporaneous element of dA. 
Furthermore, each element in dA is the beginning or the 
end of an anyon world-line. This world-line connects the 
element either to another element in dA, which has a dif¬ 
ferent time-coordinate, or to a currently existing anyon. 
Each currently existing anyon is thus connected through 
a chain consisting of strings which are alternately sub¬ 
sets of W and H to another currently existing anyon. If 
two currently existing anyons are connected this way, we 
move them towards each other, one cell per time-step, 
along the shortest possible path which is homologically 
equivalent to the chain that connects them [H] . 


The error correction problem for the fermions is much 
more involved than for the anyons because fermions are 
not only created and moved by elements of F, but also 
by elements of A and W. Fig. illustrates a process 
in which two fermions are created while correcting two 
pairs of anyons, and a process in which a fermion is 
“swallowed” during anyon error correction. Both of these 
processes illustrate that it is possible to create pairs of 
fermions which do not appear in the same measurement 
period. For this reason, the error correction problem for 
the fermions is 2 -|- 1-dimensional. We need to pair unex¬ 
pected changes of the fermionic charge which may have 
different time-coordinates. (Since before the final time- 
step we never attempt to move fermions, any change in 
fermionic charge is unexpected.) Unexpected appear¬ 
ances of fermions are due to fermion error events, or due 
to fusion of two anyons. Unexpected disappearances are 
due to an anyon “swallowing” a fermion, or due to an ap¬ 
pearance event at a location where a fermion has already 
been present. When applying MWPM to the fermionic 
problem, the weight we assign to connecting two spatio- 
temporally separated events is the 2 -|- 1-dimensional Li- 
norm. Clearly, a more sophisticated weight would take 
knowledge about anyon world-lines into account. This 
would be similar to the idea of introducing “shortcuts” in 
Ref. [57]. It may also help to weight spatial and temporal 
distances differently, and to take entropic contributions 
to the weight into account. However, we restrict to the 
Li-norm for simplicity. 

Our main result is the following theorem. 

Theorem 1. If error correction is performed as de¬ 
scribed above, there is a finite threshold Pc > 0 such that 
for p < Pc the error rate per error correction period on 
the stored quantum information is exponentially small in 

L. 

While we prove the existence of a finite lower bound 
for Pc, this lower bound is very small. Nevertheless, it 
demonstrates the important fact that the threshold is 
non-zero, and hence that continuous error correction for 
non-Abelian anyons is possible in principle. We expect 
the value for the lower bound to be very pessimistic, and 
so it should not be confused for an estimated value for 
Pc- 


III. PROOF OF THE MAIN RESULT 

Our proof is similar in nature to the proofs for the cor- 
rectability of the toric code by means of MWPM [301 HSj. 
As long as we only consider the (non-Abelian) anyons, 
the error model for them is exactly identical to the one 
for the toric code [TJ [30] with bit-flip rate p and per¬ 
fect syndrome measurements. The correctabiliy of the 
anyons thus follows from the correctability of this (very 
well-studied) error model. Indeed, Ref. [SD] contains an 
analytical proof that for this problem pc > 3.7%. 
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symbol 

natural language definition 

geometrical interpretation 

IV 

anyon world-lines 

horizontal and vertical edges of the dual lattice 

A 

anyonic (as opposed to fermionic) errors 

horizontal edges of the dual lattice 

F 

fermionic errors 

horizontal edges of the dual lattice 

dA 

unexpected changes in the anyonic charge 

cells of the lattice 

At 

anyon errors in time-step t 

horizontal edges of the dual lattice 

dAt 

unexpected changes in the anyonic charge in time-step t 

cells of the lattice 

Ht 

hypothesis about anyonic errors in time-step t 

horizontal edges of the dual lattice 

H 

all hypotheses Ht up to the present time 

horizontal edges of the dual lattice 

A’^ 

strings of anyon errors 

horizontal edges of the dual lattice 

A^ 

loops of anyon errors 

horizontal edges of the dual lattice 

At 

anyon error strings in connected component i 

horizontal edges of the dual lattice 

77. 

anyon error hypothesis in connected component i 

horizontal edges of the dual lattice 

W 

anyon world-lines in connected component i 

horizontal and vertical edges of the dual lattice 

wt 

intentional anyon movements (horizontal elements of Wt) 

horizontal edges of the dual lattice 

wr 

charge measurements detecting an anyon (vertical elements of Wt) 

vertical edges of the dual lattice 

p 

path looping around the torus 

horizontal and vertical edges of the dual lattice 

Oi 

loop formed by the disjoint union of Af and Wt 

horizontal and vertical edges of the dual lattice 

o 

union of all loops Ot 

horizontal and vertical edges of the dual lattice 

A(P) 

path P deformed by taking the symmetric difference with loop Oi 

horizontal and vertical edges of the dual lattice 

P' 

deformation of P maximizing the number of A events 

horizontal and vertical edges of the dual lattice 


FIG. 1. Natural language definition and the geometrical interpretation of all symbols that are relevant for the proof. 



FIG. 2. Two possible processes illustrating how creation and 
fusion of CT anyons can produce or “swallow” ip fermions. 
Anyonic errors are dotted, anyonic world-lines are solid, and 
fermionic world-lines are dashed. Left process: Two error 
strings produce two pairs of anyons, which are incorrectly 
paired and correspondingly brought to fusion. This process 
creates a pair of fermions with probability Right process: 
An error string creates a pair of anyons which are brought to 
fusion. Along one of the anyonic world-lines, a pair of fermions 
is created and one of the two fermions is “swallowed” by the 
nearby anyon. The second fermion is recovered only when the 
two anyons are fused. 


Our main difficulty is correcting the fermions which 
may be produced or “swallowed” during the continuous 
correction of the anyons, leading to a more involved, cor¬ 
related effective error model for the fermions. Fig. 
shows two examples of such processes. From here on, we 
consider the hypothetical completion of W. That is, we 
consider the hypothetical world-lines W that we would 
obtain if we could complete error correction in accor¬ 
dance with our hypothesis 77 at a given time and bring 


all anyons to fusion, without any further errors occurring. 
There is thus no longer a notion of “currently existing 
anyons”. Each string in W (anyon world-line) begins and 
ends at an element of dA (unexpected change in anyonic 
charge). 

The hypothetical completion of W is introduced in or¬ 
der to study whether error correction has been successful 
up to a given time, and is used to avoid explicitly model¬ 
ing a realistic fault-tolerant read-out step. We stress that 
no assumption is made that we can actually complete er¬ 
ror correction without any new errors occurring in real¬ 
ity. Similarly, we assume that the system is initially free 
of anyonic defects to avoid explicitly modeling a fault- 
tolerant read-in step. The idealized assumptions of an 
error-free initial state and a final fault-less period of er¬ 
ror correction are standard in the study of fault-tolerant 
qubit-based quantum computation (see, e.g. Ref. [45]!. 
While read-in and read-out will have to be performed in 
a fault-tolerant way in reality, we do not consider these 
for simplicity and in order to keep our results indepen¬ 
dent from the particular encoding scheme. 

We have remarked that the set A can be decomposed 
into loops and strings which connect elements in dA in 
pairs. We choose this decomposition such that each ele¬ 
ment in dA has exactly one string incident upon it. Let 
A = A® U A* be such a decomposition. (Note that the 
decomposition is in general not unique.) The sets A®, 77, 
and W can then all be decomposed into strings, each of 
which ends at an element of dA. Conversely, each ele¬ 
ment of dA has three strings incident upon it, which are 
respectively subsets of A"*, 77, and W. 















5 


Note that the sets H and W are not necessarily dis¬ 
joint: it can happen that we attempt to move an anyon 
to an adjacent cell and, between the same two rounds of 
charge measurement, an anyonic error affecting the same 
two cells happens. Let us thus study the disjoint union 

A^UHUW = 

{(e, a) : e £ A®} U {(e, /i) : e £ H} U {{e,w) : e £ W} . 

( 2 ) 


Furthermore, each error event can create at most two 
anyons, so 

h<\At\. ( 6 ) 

Recall from Eq. (|^ that 

(7) 

Combining the above inequalities, we find 


Here, the index i £ {a,h,w} in the ordered pair (e,i) 
tells us which of the three sets H®, H, or W the edge e 
belongs to. If, for example, an edge is an element of both 
H and W, there will thus be two corresponding edges in 
A'^UHUW. 

The set A^ Li H LiW forms a trivalent graph, with each 
vertex corresponding to an element of dA, and having 
an H®, an H-, and a IT-string incident upon it. Let us 
study minimal connected components of this graph. Let 
A®, Hi, and Wt denote the union of all strings in H®, H, 
and W, respectively, that belong to connected component 
i. Finally, let Wi denote the set of horizontal edges (i.e., 
intentional anyon movements) and the set of vertical 
edges (i.e., charge measurements which detect an anyon) 
in Wi- Recall that Af and Hi consist of horizontal edges 
only. 

From the way our error correction procedure is defined, 
we have the inequalities 

\wt\< m < . ( 3 ) 

The first inequality is due to us moving anyons along the 
shortest path which is homologically equivalent with Hi. 
We could always choose \ Wl^\ = \Hi\ by undoing exactly 
the errors that happened according to our hypothesis. 
The second inequality is due to using MWPM for error 
correction. Assume by contradiction that \Hi\ > |A®|. 
Then, replacing Hi with Af would yield a perfect match¬ 
ing of the unexpected changes in anyonic charge which 
is of lower weight than the one returned by the MWPM 
algorithm, which contradicts its definition. 

Now consider the loops Oi = Af Li Wi. The following 
lemma is the main technical tool that we use in order to 
deal with these. 

Lemma 1. All loops Oi satisfy 

|0.|<4|A®|. (4) 

Note that for the simplest possible process, a single 
anyon error event that is immediately corrected (|A®| = 
1, m = 1, \Wf\ = 1, IWfl = 2, |0,| = 4), the bound is 
tight. 

Proof. Let fi denote the number of fusion events of a 
pair of anyons in Oi. An anyon needs to be moved away 
from each location at which it appears. However, before 
fusion two anyons may be adjacent so that we need to 
move only one of them. We thus have 

\wt\>\wn-f^. ( 5 ) 


|W;|<2|A®|. (8) 

For the total length of the loop, we find, using Eqs. Q 
and (|^, 

|0.| = |A®| + |H^'*| + |Wr|<4|A®|. (9) 

□ 


The following lemma provides a necessary condition 
for the failure of error correction. 


Lemma 2. A failure of error correction requires a ho¬ 
mologically non-trivial closed path P (a loop) satisfying 

7\PnA\ + \PnF\>\p\/2. (10) 


Proof. Let us first study the possibilities for error correc¬ 
tion failing for the anyons (as opposed to the fermions). 
Recall that we have decomposed the set of anyonic errors 
A into loops A’’ and strings A®. If one of the loops which 
are subsets of A^ is homologically non-trivial, Eq. (10) 


will obviously be satisfied, as we can choose P to be the 
corresponding loop and have |P n A| = |P|. The second 
possibility for error correction for the anyons failing is 
that one of the loops Oi is topologically non-trivial. In 
this case, we choose P — Oi and are done, since by use 
of Lemma [T] we have 


7\0,nA\ = 7\Af\>l\0,\. 


( 11 ) 


So assume from now on that all loops which are subsets 
of A* and all loops Oi are homologically trivial, and that 
error correction failing is due to the fermionic part of the 
problem. 

Clearly, MWPM failing to correct the fermions requires 
that there be a homologically non-trivial closed path P 
containing at least |P|/2 edges that have been affected 
by an event that can possibly have created or moved 
fermions, for otherwise the minimum-weight correction 
of the fermions will never move a fermion around the 
torus. We assume pessimistically that each edge in A* 
and in O = (J - Oi (anyon error event or anyon world-line) 
counts as a potential fermion error event. So formally, we 
need a path P with 

|Pn (A'uOUF)I > |P|/2. (12) 


We will prove that if there is such a path P, there is a 
(possibly identical) path P' which is homologically equiv¬ 
alent to P and satisfies the inequality in the lemma, i.e., 
7|P'n A| -b |P'nP| > |P'|/2. 
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Given a loop Oi with n P ^ 0, we can consider 
the “deformed” path Di{P) = (P \ Oi) U (Oi \ P). The 
path P' is obtained by applying a (possibly empty) set 
of deformation operations Di to P. Since all of the loops 
Oi are homologically trivial, the deformed path P' will 
be homologically equivalent to P. We define the path P' 
such that the number of A events in the path is maxi¬ 
mized; i.e., such that 

\p' n = max{|p n |A(P) n , (13) 

for all loops with Oi n P 7 ^ 0. Equivalently, the path P' 
is defined such that 


IAn^^'l<l^-nP'|. (14) 

Let us define = A^\0 and F = F\0. By assump¬ 
tion, we have 

0 < |Pn (A' uouP)| - |P|/2 
= |p n i'l -b |P n 0| -b |P n P| - |p|/2 
= iPni'l -b |PnO| -b |PnP| 
-(|Pn0|/2+|P\0|/2). (15) 


Note that P\0 is not affected by deformation operations, 
i.e. P\0 — P' \ 0. Since P n P C P \ O, we also have 
PnF = P' n P and similarly PH A'' = P'HAK Therefore 


0 < |p' n i'l + |p n o|/2 + |p' n p| - |p' \ o|/2 
= iP'ni'l + (|Pn 0| + |P'n 0|)/2 + |p'np| 
-(|P'\0|/2+|P'n0|/2) 

= |P'ni'| + ^(|Pno,| + |P'no,|)/2 

i 

+ |p'np| - |p'|/2. (16) 


If P n Oi = P' n Oi we find 


(|pno,| + |p'no,|)/2 
= |Pno,| 

= \PnAt\ + \Pno,\At\. (17) 

Since P n Oi \ blf C Oi \ A| and Af C Oi we have 


\pno,\At\<\o,\At\ = \o,\-\At 


(18) 


Combining Eqs. ( [T7| ) and (18) with Lemma we arrive 
at 


(|Pno,| + |P'no,|)/2< |Pnblf| + 3|blf|. (19) 


Using Eq. (14), we obtain 

= lA* n P'l + lA,* \ P'l < 2 |bl^ n P'l = 2 |blf n p| . 

( 20 ) 

We finally find 

(|P n 0,1 + \P' n o ,|)/2 < 7|P' nA^\. (21) 


If, on the other hand, P H Oi ^ P' H Oi, we find, using 
Lemma for the first inequality, 

(|Pn 0,1 + |P'n 0|,)/2 = |0,|/2<2|bl^| <4|P'U bill. 

( 22 ) 


So in both cases Eq. (21) holds and we find from Eq. (16) 
that 


0 < |P' n i'l -b 7^ |p' n All + |P' n P| - |p'|/2 

i 

< 7|p'n A| + |P'nP| - |P'|/2. (23) 


□ 

Theorem directly follows from the following lemma. 

Lemma 3. The probability per time-step of a path as 
in Lemma is exponentially suppressed with L if p < 
15-14 « 3 X 10-11'. 


Proof. Consider two lines of length L looping in homo¬ 
logically non-equivalent ways around the torus. Path P 
needs to cross at least one of them. Since the two lines 
can be crossed at 0{L) locations, and a path in a three- 
dimensional cubic lattice can at each step turn into 5 
directions, there are at most 5 ^“'- 0 (iog(L)) dosed paths of 
length i in the lattice crossing any of the two lines at a 
given time. Let n = |P n A| + |P n P| be the number of 
error events along the path, and let £ = |P|. For a path 
satisfying 7|P n A| + |P n P| > |P|/2, we need 14n > £. 
With fixed locations of the n errors, the probability of 
such a path is at most p” < p^^i^. In a path of length 
i, there are no more than 3^ possibilities for picking the 
locations of A and P events. The probability per time- 
step of a path satisfying 7|P n A| + |P n P| > |P|/2 is 
thus upper-bounded by 

OO 

5 ^+o(iog(L)) 3 y/i 4 ^ (24) 

t=L 

which is exponentially suppressed with L if < 

1 . □ 


IV. CONCLUSIONS 

Topological quantum computing holds the promise of 
processing quantum information with hardware which 
has intrinsically low error rates. Still, these rare errors 
need to be corrected in a truly large-scale computation. 
Here, we have demonstrated the feasibility of this task 
for Ising anyons under the realistic assumption that er¬ 
rors keep happening as we correct those from previous 
rounds. 

Note that another approach towards error correction 
would be to use qubits encoded in non-Abelian anyons 
as physical qubits for a further round of error correc¬ 
tion. This would lead to a standard decoding problem. 
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although perhaps with a rather complex error model due 
to the underlying topological processes. In such a case, 
the available gate set on the final logical qubits would 
be restricted by the code used in the additional round 
of error correction. This would remove some of the ad¬ 
vantages of using non-Abelian anyons, and may carry a 
large resource overhead in comparison to standard pro¬ 
posals for physical qubits. 

One would hope for a threshold proof for further non- 
Abelian anyon models, including those for which MWPM 
cannot be applied to perform error correction [55H1H]. 
Unfortunately, the highly general proof in Ref. [53] does 
not allow for straightforward generalization to the contin¬ 
uous case. Furthermore, it is an open problem to study 
error correction for non-Abelian anyons with charge mea¬ 
surements that can give an incorrect result. 

Finally, it would be valuable to get a better idea of the 
“true” threshold for our setup. Given the crudeness of 
our arguments, we expect our threshold of Pc ~ 3 x 10“^^ 


to be rather pessimistic. A better estimate of the true 
threshold value Pc could be obtained via a more inge¬ 
nious analytical approach, or by numerical simulations, 
extending the work of Ref. [5^ to the continuous case. 
The “true” thresholds for Abelian models can often be 
assessed by finding the phase-transition in a related clas¬ 
sical statistical mechanics model [SOI SZl SH]. Whether 
something similar can be done for non-Abelian models 
remains an open problem. 
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